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We investigate the properties of the hadron-quark mixed phase in compact stars using 
a Brueckner-Hartree-Fock framework for hadronic matter and the MIT bag model for 
^ ' quark matter. We find that the equation of state of the mixed phase is similar to that 

' given by the Maxwell construction. The composition of the mixed phase, however, is very 

^ different from that of the Maxwell construction; in particular, hyperons are completely 

, suppressed. 
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. 1. Introduction 

o 

Matter in neutron stars (NS) has a variety of density and chemical component due 
. ^ , to the presence of gravity. At the crust of neutron stars, there exists a region where 

' the density is lower than the normal nuclear density, po « 0.17 fm""^ over a couple of 

^ . hundreds meters. The pressure of such matter is mainly contributed by degenerate 

electrons, while baryons are clusterized and have little contribution. Due to the 
gravity pressure and density increase in the inner region (in fact, the density at 
the center amounts to several times po). Cold catalyzed matter consists of neutrons 
and the equal number of protons and electrons under chemical equilibrium. Since 
the kinetic energy of degenerate electrons is much higher than that of baryons, 
the electron fraction (or the proton one) decreases with increase of density and 
thus neutrons become the main component and drip out of the clusters. In this 
way baryons come to contribute to the pressure as well as electrons. At a certain 
density, other components such as hyperons and strange mesons may emerge. At 
even higher density, hadron-quark deconfinement transition may occur and quarks 
in hadrons are liberated. 

It is well known that hyperons appear at several times po and lead to a strong 
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softening of the EOS with a consequent substantial reduction of the maximum 
neutron star mass. Actually the microscopic Brueckner-Hartree-Fock approach gives 
much lower masses than current observation values of ^ IAMq. 

On the other hand, the hadron-quark deconfinemcnt transition is believed to 
occur in hot and/or high-density matter. Taking EOS of quark matter within the 
MIT bag model, the maximum mass can increase to the Chandrasekhar limit once 
the deconfinement transition occurs in hypcron matter. -'^■^ 

The deconfinemcnt transition from hadron to quark phase may occur as a first- 
order phase transition. There, a hadron-quark mixed phase should appear, where 
charge density as well as baryon number density is no more uniform. Owing to the 
interplay of the Coulomb interaction and the surface tension, the mixed phase can 
have exotic shapes called pasta structures.^ 

Generally, the appearance of mixed phase in matter results in a softening of the 
EOS. The bulk Gibbs calculation (BG) of the mixed phase, without the effects of 
the Coulomb interaction and surface tension, leads to a broad density region of the 
mixed phase (MP).^ However, if one takes into account the geometrical structures 
in the mixed phase and applies the Gibbs conditions, one may find that MP is con- 
siderably limited and thereby the EOS approaches to the one given by the Maxwell 
construction (MC).^'^''^ 

In this report we explore the EOS and the structure of the mixed phase dur- 
ing the hyperon-quark transition, properly taking account of the Gibbs conditions 
together with the pasta structures. 



2. Numerical Calculation 

The numerical procedure to determine the EOS and the geometrical structure of 
the MP is similar to that explained in detail in Ref. 3. We employ a Wigner-Seitz 
approximation in which the whole space is divided into equivalent Wigner-Seitz cells 
with a given geometrical symmetry, sphere for three dimension (3D), cylinder for 
2D, and slab for ID. A lump portion made of one phase is embedded in the other 
phase and thus the quark and hadron phases are spatially separated in each cell. A 
sharp boundary is assumed between the two phases and the surface energy is taken 
into account in terms of a surface-tension parameter a. The energy density of the 
mixed phase is thus written as 



1 

Vw 



(1) 

where the volume of the Wigner-Seitz cell Vw is the sum of those of hadron and 
quark phases Vh and Vq, S the quark- hadron interface area, en, eg and are en- 
ergy densities of hadrons, quarks and electrons, which are functions of local densities 
/9a(r) (a = n,p,A,T,~,u,d,s,e). The Coulomb potential Vq is obtained by solving 
the Poisson equation. For a given density ps, the optimum dimensionality of the 
cell, the cell size Rw, the lump size R, and the density profile of each component 
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are searched for to give the minimum energy density. We employ cr = 40 MeV/ fm^ 
in the present study. 

To calculate en in the hadron phase, we use the Thomas-Fermi approximation 
for the kinetic energy density. The potential-energy density is calculated by the 
nonrelativistic BHF approach^ based on microscopic NN and NY potentials. 



2— n,p,A,S 
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(3) 
(4) 

(5) 



The interaction parameters are chosen to reproduce the scattering phase shifts. 
Nucleonic three-body forces are included in order to (slightly) shift the saturation 
point of purely nucleonic matter to the empirical value. 

For the quark phase, we use the MIT bag model with massless u and d quarks and 
massive s quark with = 150 MeV. The energy density eg consists of the kinetic 
term by the Thomas-Fermi approximation, the leading-order one-gluon-exchange 
term^ proportional to the QCD fine structure constant as, and the bag constant 
B. Here we use B — 100 MeV/fm'^ and = to get the quark EOS which crosses 
the hadronic one at an appropriate baryon density. 



3. Hadron-Quark Mixed Phase 
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Fig. 1. Left: Density profiles and Coulomb potential Vc within a 3D (quark droplet) Wigner- 
Seitz cell of the MP at pg = 0.4 fm~^. The cell radius and the droplet radius are Ryy = 26.7 
fm and R = 17.3 fm, respectively. Center: Same as the left panel for MC case. The radius r is in 
arbitrary unit since there is no specific size. Right: The case of BG calculation. 



Figure [T] illustrates an example of the density profile in a 3D cell. One can 
see the non-uniform density distribution of each particle species together with the 
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finite Coulomb potential; charged particle distributions are rearranged to screen the 
Coulomb potential. The quark phase is negatively charged, so that d and s quarks 
are repelled to the phase boundary, while u quarks gather at the center. The protons 
in the hadron phase are attracted by the negatively charged quark phase, while the 
electrons are localized to the hadron phase. This density rearrangement of charged 
particles causes the screening of the Coulomb interaction between phases. 

In the center and right panels, compared are the cases of MC and BG. MC 
assumes the local charge neutrality, while the BG does not. One can see that the 
local charge neutrality in the full calculation lies between two cases. The localization 
of electrons which is one of the charge screening effects, reduces the local charge 
density. But the local charge density remains still finite to some extent. 




0.2 0.4 , 0.6 

PB[fm-3] 



Fig. 2. EOS of the MP (thick curves) in comparison with pure hadron and quark phases (thin 
curves), hadron (pg < 0.44 fm""^) or quark (pg > 0.44 fm""^) phases. Each segment of the MP is 
chosen by minimizing the energy. 

Figure [2] compares the resulting EOS with that of the pure hadron and quark 
phases. The thick black curve indicates the case of the MC, while the colored lines 
indicate the MP in its various geometric realizations starting with a quark droplet 
structure and terminating with a bubble structure. Note that the charge screening 
effect always tends to make matter locally charge neutral to save the Coulomb 
interaction energy. Hence, combined with the surface tension, it makes the non- 
uniform structures mechanically less stable and limits the density region of the 
MP."^ Consequently the energy of the MP is only slightly lower than that of the 
MC. 

However, the structure and the composition of the MP are very different from 
those of the MC, which is demonstrated in Fig. [3l where we compare the particle 
fractions as a function of baryon density in the full calculation (left panel) and the 
MC (right panel). One can see that the compositions are very different in two cases. 
In particular, a relevant hyperon (S^) fraction is only present in the MC. 

The suppression of hyperon mixture in the MP is due to the absence of the 
charge- neutrality condition in each phase. As shown in Fig. [H hyperons (E~) ap- 
pear in charge-neutral hadronic matter at low density (0.34 fm~^) to reduce the 
Fermi energies of electron and neutron. Without the charge-neutrality condition. 
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Fig. 3. Particle fractions in the MP by the full calculation (left panel) and the MC (right panel). 




Fig. 4. Upper panel: Particle fractions of neutral matter with electrons (corresponding to neutron 
star matter). Lower panel: The same quantity for charged matter without electrons, the low-density 
part of which corresponds to symmetric nuclear matter. Both cases require chemical-equilibrium 
condition. 



on the other hand, there appears symmetric nuclear matter at lower density and 
hyperons will be mixed above 1.15 fm""^ due to the large hyperon masses. The MP 
has positively charged hadron phase and negative quark phase though the Coulomb 
screening effect diminishes the local charge density. This brings the feature of sym- 
metric nuclear matter into the hadron phase and, consequently, the mixture of 
hyperons is suppressed.^ 



4. Neutron Star Structure 

Having the EOS comprising hadronic, mixed, and quark phase in the form P{e), the 
equilibrium configurations of static NS are obtained in the standard way by solving 
the Tolman-Oppenheimer-Volkoff (TOV) equations^ for the pressure P{r) and the 
enclosed mass m(r), 

dP_ _ Gme (1 + P/e) (l + Anr^P/m) 
^ ^ 1 - 2Gm/r ' 

— = Anr^e , 7) 
dr 

being G the gravitational constant. Starting with a central mass density e(r = 0) = 
Cc, one integrates out until the surface density equals the one of iron. This gives the 
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Fig. 5. (Color online) Neutron star mass-radius relations for different EOS and three different 
hadron-quark phase transition constructions. For the hybrid stars (blue and black curves), the 
dashed lines indicate the Maxwell (upper curves) or bulk Gibbs (lower curves) constructions and 
the solid lines the mixed phase of the full calculation. 




Fig. 6. (Color online) fnternal structure of a 1.4 Mq neutron star obtained with three different 
phase transition constructions. The upper panels show total energy density and pressure and 
the lower panels the overall particle fractions as functions of the radial coordinate of the star, 
using the bulk Gibbs calculation (left panel), the mixed phase of the full calculation with cr = 
40 MeV/fm^ (central panel), and the Maxwell construction (right panel). In all cases cts = and 
B = 100 MeV/fm-3 are used. 



stellar radius R and its gravitational mass M = m(R). For the description of the NS 
crust, we have joined the hadronic EOS with the ones by Negele and Vautherin^" 
in the medium-density regime, and the ones by Feynman-Metropolis-Teller^^ and 
Baym-Pethick-Sutherland^^ for the outer crust. 

Fig. E] compares the mass-radius relations obtained with the different models. 
The purely nucleonic EOS (green curve) yields a maximum NS mass of about 
1.82 Mq, which is reduced to 1.32 Mq when allowing for the presence of hyperons 
(red curve). This feature has been shown to be fairly independent of the nucleonic 
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and hyperonic EOS that are used.^'^ The canonical NS with mass of about 1.4 
can therefore not be purely hadronic stars in our approach. In fact, the inclusion of 
quark matter augments the maximum mass of hybrid stars to about 1.5 A/q: 

In general, the Maxwell construction leads to a kink in the M{R) relation, 
because the transition from a hadronic to a hybrid star occurs suddenly, involving a 
discontinuous increase of the central density when the quark phase onsets in the core 
of a star. The bulk Gibbs calculation yields smooth mass-radius relations involving 
a continuous transition from a hadronic to a hybrid star beginning at rather low 
central density corresponding to very low NS mass. 

The MP construction by the full calculation lies between the two extreme cases, 
and with our choice of cr = 40 MeV/fm^ it is rather close to the Maxwell construc- 
tion, smoothing out the kink of the hadron-hybrid star transition. This transition 
occurs generally at a fairly low NS mass, even below the natural minimum mass 
limit due to the formation via a protoneutron star^** and is thus an unobservable 
feature. 

Whereas the maximum masses are practically independent of the phase transi- 
tion construction, there are evidently large differences for the internal composition 
of the star. This is illustrated in Fig.[5]which show the total energy density, pressure, 
and particle fractions as a function of the radial coordinate for a 1.4 Mq NS. One 
observes with the bulk Gibbs construction (left panels) a coexistence of hadrons 
and quarks in a significant range of the star, whereas with the MC (right panels) an 
abrupt transition involving a discontinuous jump of energy and baryon density oc- 
curs at a distance r « 7.5 km from the center of the star. The small contamination 
with S~ hyperons in the hadronic phase is not visible on the scale chosen. The MP 
with the full calculation (central panels) lie between the two extreme cases, hadrons 
and quarks coexisting in a smaller range than in the bulk Gibbs cases. 

5. Summary 

In this article we have studied the properties of the mixed phase in the quark 
deconfinement transition in hyperonic matter, and their influence on compact star 
structure. The hyperonic EOS given by the BHF approach with realistic hadronic 
interactions is so soft that the transition density becomes very low if one uses the 
MIT bag model for the quark EOS. 

The hyperon-quark mixed phase was consistently treated with the basic thermo- 
dynamical requirement due to the Gibbs conditions. We have seen that the resultant 
EOS is little different from the one given by the Maxwell construction. This is be- 
cause the finite-size effects, the surface tension, and the Coulomb interaction tend 
to diminish the available density region through the mechanical instability, as has 
also been suggested in previous articles. 

For the bulk properties of compact stars, such as mass or radius, our EOS gives 
similar results as those given by the Maxwell construction. The maximum mass 
of a hybrid star is around 1.5 Mq, larger than that of the purely hyperonic star. 
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« 1.3 Mq. Hence we may conclude that a hybrid star is still consistent with the 
canonical NS mass of 1.4 Mq, while the masses of purely hyperonic stars lie below 
it. 

On the other hand, the internal structure of the mixed phase is very different; 
e.g., the charge density as well as the baryon number density are nonuniform in the 
mixed phase. We have also seen that the hyperon number fraction is suppressed in 
the mixed phase due to the relaxation of the charge-neutrality condition, while it 
is always finite in the Maxwell construction. This has important consequences for 
the elementary processes inside compact stars. For example, coherent scattering of 
neutrinos off lumps in the mixed phase may enhance the neutrino opacity. Also, 
the absence of hyperons prevents a fast cooling mechanism by way of the hyperon 
Urea processes.^* These results directly modify the thermal evolution of compact 
stars. 
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